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摘 要 : 为 揭示 不 同 土地 类 型 蒸 散 发 量 ( ET) 差异 特征 ,基于 SEBS (surface energy balance system) 模 型 ,结合 气象 数 
据 . 遥 感 影像 及 土地 利用 分 类 结果 ,对 河套 灌区 2017 年 生长 季 日 57 进行 SEBS 模型 估算 ,分 析 灌 区 ET 时 空 变化 特 


征 及 不 同 土地 类 型 的 7 差异。 结果 表明 :Q@ 通过 实测 数据 和 FAO Penman-Monteith( PM ) 计算 值 对 SFEBS 模型 估算 


结果 精度 验证 , 蒸 渗 仪 观测 数据 、 涡 度 相 关系 统 (EC) 波 文 比 能 量 平衡 系统 (BREB) 和 PM 计算 值 的 相对 误差 平均 
值 分 别 为 :16.66% 、18.34% 、14.93% 、14.03% ,表明 SEBS 模型 可 以 用 来 估算 河套 灌区 ET, O 生长 季 内 ,灌区 日 


ET 大 小 依次 为 :7 月 > 6 月 >5 月 >8 月 >9 
日 £7 最 小 ,为 3.51 mm: d7 


Jo A ET 呈 单 峰 形变 化 ,7 H H ET 达到 最 大 值 5.21 mm .4d”,9 H 
'。@@ 不 同 土地 类 型 的 日 £7 大 小 依次 为 :水 体 > 林 地 > 耕地 > 草地 > 城市 > 荒地 。 灌 
区 内 ET 分 布 与 土地 类 型 具有 高 度 一 致 性 ,总 体 分 类 精度 和 Kappa 系数 分 别 为 80. 63% F1 0.75, O 在 生长 季 的 不 
同时 间 ,水 体 ET 均 处 于 较 高 值 ,7 H H ET 达到 最 大 值 6.034 mm + d ,9 


a 


日 £7 值 最 小 ,为 4.177 mm: d7', 其 


次 林地 日 67 较 大 ,7 月 日 £7 值 为 5.977 mm + d ;9 月 日 ET7 值 为 4.147 mm d`, BEA ET 变化 趋势 与 作物 生 
长 情况 相 一 致 ,作物 生长 旺盛 期 ET (AGA 5.851 mm . d ,作物 收获 期 ET 值 为 3.952 mm - d`! 
关键 词 : HRA; AMA; SEBS 模型 ; Softmax 分 类 方法 ; 土地 类 型 ; 河套 灌区 


蒸 散发 (ET7) 是 地 表 水 量 和 热量 平衡 的 重要 参 
E ,是 水 分 吸收 能 量 从 液态 转变 为 气态 的 过 程 ,包含 


结果 显示 SEBS 模型 的 结果 要 优 于 蒸 散 发 产品 。 王 
瑶 ” 利用 SEBS 模型 分 析 黑 河中 游 作物 需 水 量 的 时 


土壤 水 体 和 植被 的 蒸发 .蒸腾 ,也 是 衡量 植物 生长 
状况 和 作物 产量 的 重要 指标 o RPR æki 
环 和 能 量 循环 的 重要 组 成 部 分 ,影响 着 区 域 乃至 全 
球 的 生态 .气候 及 农业 生产 ,所 以 对 藻 散 发 进行 有 效 
的 估算 是 十 分 必要 的 。E7 估算 方法 研究 已 有 200 
多 年 的 历史 ,传统 的 ET 估算 是 基于 点 位 的 地 面 观 
测 , 只 适用 于 林地 、 农 地 等 小 尺度 范围 。 考 虑 到 大 区 
人 
将 点 位 观测 结果 推广 到 区 域 尺度 。 近 年 来 ,遥感 
技术 具有 大 范围 实时 高效、 费用 低廉 等 优点 ,为 估 


算 大 区 域 的 地 表 ET 开辟 了 一 条 新 途径 3] 。 遥 感 技 
术 可 以 在 时 间 和 空间 上 快速 获取 大 面积 的 地 物 光 谱 


信息 ,能 及 时 反映 出 由 于 水 分 盘 亏 所 引起 的 地 表 ET 
变化 。 王 卫 光 等 外 使 用 SEBS 模型 结合 MODIS 数 
据 与 现 有 的 MOD16 .ETzhang 燕 散发 产品 进行 对 比 ， 


收 稿 日 期 : 2019-09-04; ”修订 日 期 : 2019 -12 -02 
基金 项 目 : 国家 自然 科学 基金 项目 (51769021 ,51839006 ) 


© 


作者 简介 : EIE (1993 — ) , 女 ,硕士 研究 生 , 研 究 方向 为 农业 水 利和 遥感 信息 


通讯 作者 : 李 瑞 平 


.了 -mail:nmghp@ 163. com 


空 分 布 格局 ,并 对 SEBS 模型 的 输入 参数 进行 了 敏 
感性 分 析 。 张 圆 等 中 采用 SEBS 模型 对 新 疆 昌 吉州 
呼 图 壁 县 蒸 散发 量 进 行 估 算 , 反 映 了 天 山北 坡 典型 
县 域 蒸 散 量 的 时 空 变化 特征 。 伴 随和 遥感 数据 大 量 增 
多 , 且 精 度 逐 步 提高 ,针对 河套 灌区 ET 研究 有 了 更 
广泛 的 素材 。 张 稚 等 “2009 年 基于 MODIS 数据 通 
过 模型 反 演 河套 平原 2006 年 日 87, 并 对 河套 平原 
ET 进行 了 分 析 。 杨 十 亭 等 "2012 年 利用 双 源 模型 
对 河套 平原 潜在 ET 进行 了 对 比 研究 ,通过 假设 不 
同 的 植被 覆盖 状况 ,较为 详细 地 比较 了 3 PRIORE 
散发 模型 在 估算 潜在 ET 以 及 区 分 潜在 蒸发 量 与 淤 
在 腾 发 量 时 的 模拟 效果 。 这 些 研 究 为 估算 区 域 ET 
提供 了 可 靠 适 用 的 理论 基础 。 

内 蒙古 河套 灌区 属 干 旱 气候 区 域 , 水 资源 短缺 ， 
Mt E BES , 且 区 域 分 布 不 均衡 ,农业 生产 主要 靠 引 
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黄河 水 灌溉 。 目 前 ,虽然 针对 河套 灌区 ET 研究 已 
有 部 分 成 果 ,但 难以 准确 反映 出 整个 灌区 ET 空间 
异 质 性 。 土 地 利用 的 方式 对 灌区 ET 也 有 显著 的 影 
i, ET 随 不 同 土地 类 型 而 异 。 明 确 ET 实际 时 空格 
局 变化 特征 ,对 于 合理 分 配水 资源 .提高 水 资源 利用 
率 具有 重要 意义 。 本 文 基于 SEBS 模型 估算 灌区 
E7, 定 量 研 究 不 同 土地 类 型 生长 季 的 ET 时 空 变化 
规律 及 其 差异 特征 ,明确 认识 土地 利用 对 E7 的 影 
响 。 对 于 调节 灌区 气候 保护 生态 环境 ,促进 农业 生 
产 有 着 重要 的 现实 意义 ,为 进一步 改善 灌区 土地 合 
理 利用 、 有 效 分 配水 资源 提供 一 定 的 技术 支持 。 


1 研究 区 概况 及 数据 来 源 


1.1 研究 区 概况 

Py Ser vel eK a FE ER REO a 
后 旗 .临河 区 AL .乌拉 特 前 旗 5 个 地 区 ,东西 约 
250 km ,南北 约 50 km。 北 依 阴 山 山 脉 的 狼 山 .乌拉 
山南 葛 洪 积 扇 , 南 临 黄河 , 东 至 包头 市 郊 , 西 接 乌 兰 
布 和 沙漠 , 属 中 温带 季风 气候 ,地 理 坐 标 40°12" ~ 
41°18'N,106°20' ~ 109°29'E。 灌 区 地 势 平坦 ,西南 
高 ,东北 低 ,海拔 高 度 1 007 ~ 1 050 m, 坡 度 
0.125%o ~0.2 %o。 实际 灌溉 面积 约 6. 01 x 10° 
hm?!’ , 占 总 土地 面积 的 51.2% 左右 , 引 黄 水 量 占 灌 
区 总 耗 水 量 的 86. 6% ,是 灌区 最 重要 的 水 分 来 源 。 
灌区 降水 量 少 , 荧 发 量 大 ,年 平均 降水 量 140 ~ 230 
mm ,年 平均 蒸发 量 1 900 ~ 2 300 mm。 夏 季 高 温 干 
旱 冬季 严寒 少 雪 ,无 霜 期 短 ,封冻 期 长 ,全 年 以 西风 
和 西北 风 为 主 , 日 温差 大 ,日 照 时 间 长 。 研 究 区 高 程 
及 空间 位 置 如 图 1 所 示 。 
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1.2 数据 来 源 
1.2.1 BR 遥感 影像 是 SEBS 模型 进行 区 
SET 估算 的 必要 条 件 ,Landsat 8 遥感 影像 从 地 理 
28 [a] ACHE z (http://www. gscloud. en/) 下载。 利用 
EZR Be SEBS 模型 反 演 区 域 57, 对 影像 数据 的 
质量 要 求 较 高 ,综合 考虑 研究 区 作物 种 植 时 间 及 生 
长 情况 ,选取 河套 灌区 2017 年 生长 季 5 一 9 月 晴天 
条 件 且 云 量 较 少 .图 像 质 量 较 好 的 遥感 影像 进行 ET 
分 析 。 对 影像 数据 分 别 进行 辐射 定 标 、FLASH 大 气 
校正 等 预 处 理 , 反 演出 SEBS 模型 ,计算 ET 所 需 的 
遥感 参数 ,包括 地 表 温 度 、 地 表 比 辐射 率 、 地 表 反 照 
率 归 一 化 植被 指数 植被 覆盖 度 . 叶 面 积 指数 、 植 被 
高 度 .位移 高 度 .表面 粗糙 度 。 
1.2.2 气象 数据 及 其 他 辅助 数据 地 面 气象 数据 
来 源 于 中 国 气 象 数 据 网 (http ://data. cma. en/) 及 部 
分 当地 气象 站 。 收 集 处 理 2017 年 研究 区 内 及 周边 
的 14 个 气象 站 5 一 9 月 日 值 气象 数据 ,气象 站 分 布 
如 图 1 所 示 。 其 中 位 于 研究 区 域外 的 9 个 气象 站 主 
要 用 于 局 部 空间 插值 ,以 保证 研究 区 边界 气象 参数 
的 空间 插值 结果 连续 平滑 。 气 象 数 据 包 括 日 照 时 
数 \ 日 平均 气压 、 日 地 表 压 力 、 日 平均 气温 、 相 对 湿 
度 .日 平均 风速 、 比 湿 。 为 提高 SEBS 模型 估算 精 
度 ,将 不 同 气象 参数 进行 Kriging 空间 插值 ,得 到 与 
遥感 影像 空间 分 辨 率 相同 的 气象 栅 格 数据 。 

精度 验证 数据 分 为 实测 数据 和 计算 数据 两 类 。 
以 FAO Penman-Monteith(PM) 计 算 的 灌区 内 5 个 气 
象 站 点 的 ET 值 作 为 理论 数据 ,实测 数据 源 自 2017 
年 5 一 9 月 临河 试验 地 的 蒸 渗 仪 观测 ET 数据 和 分 
子 地 涡 度 相关 (EC ) 及 波 文 比 能 量 平 衡 (BREB ) A 
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Fig. 1 


图 1 研究 区 高 程 及 空间 位 置 示意 图 


Schematic diagram of research area elevation and location 
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THR 
统 得 到 的 ET Eo 
(1) 蒸 渗 仪 的 观测 数据 获取 
临河 试验 地 实测 ET 值 采 用 微型 蒸 渗 仪 测定 ， 


测定 设置 因素 分 别 为 材料 和 尺寸 。 其 中 ,内 外 桶 材 
料 均 采用 PVC 材质 ,内 桶 直径 10 cm, 外 桶 直径 15 
cm ,内 外 桶 深度 均 为 20 cm, 封 底 材料 选用 细 纱 布 。 
于 试验 期 间 每 天 08 :00 采用 精度 0.1 g 的 电子 天 平 
(PL 6001 - 工 ,Mettler Toledo Inc. ,瑞士 ) 对 微型 蒸 涂 
仪 称 重 ,根据 两 次 重量 差 ,经 面积 换算 再 加 上 灌溉 或 
降水 量 得 到 每 天 的 ET 值 。 为 保证 测量 数据 的 准确 
性 ,在 田间 不 同位 置 布设 7 个 相同 规格 的 微型 蒸 渗 
仪 , 取 其 平均 值 。 

(2) 涡 度 相关 及 波 文 比 能 量 平衡 系统 的 观测 数 
据 获 取 

临河 分 子 地 安装 了 涡 度 相关 及 波 文 比 能 量 平衡 
系统 ,于 2017 年 5 一 9 月 底 连 续 测量 , 均 位 于 被 同一 
作物 所 包围 的 观测 场 ,为 系统 测量 ET 值 提供 了 足 
够 的 获取 长 度 。 以 自立 式 铁塔 作为 能 量 观测 塔 ,分 
别 架设 涡 度 相关 及 波 文 比 能 量 平衡 系统 ,在 作物 最 
大 冠 层 高 度 上 方 1.5 m 处 安装 带 有 开放 式 红外 H,O 
和 CO, 分 析 仪 的 CSAT3 三 维 超声 波 风速 计 (LI - 


ChinaXiv 合 作 期 刊 


7500 ,LI-COR Inc. ,美国 ) , 测 水 平和 垂直 方向 风速 
的 脉动 值 ,并 在 作物 生育 期 经 过 一 次 位 置 调 整 。 
TCAV 热电 侦探 头 测 土 壤 的 平均 温度 ,CS616 土壤 
湿度 传感器 测 土壤 水 分 。 在 观测 塔 底部 装 有 观测 小 
屋 ,用 于 观测 数据 的 读 取 .记录 及 临时 处 理 。 通 量 塔 
一 侧 的 传感器 通过 线路 与 观测 小 屋 相 连接 ,数据 以 
10 Hz 采样 并 使 用 CR3000 数据 记录 仪 以 30 min 间 
隔 记录 。 在 通 量 塔 的 另 一 侧 , 架 设 高 度 为 最 大 树冠 
上 方 1 m 处 , 装 有 CNR4 - L10 44h Bit (Kipp & 
Zonen ,荷兰 ) 用 于 测量 向 上 和 向 下 的 长 短波 净 辐 射 
及 日 照 时 数 。 在 土壤 表面 以 下 约 6 cm 处 埋设 两 个 
HFP01 热流 板 (Hukseflux ,荷兰 ) , 测 一 段 时 间 内 的 
地 热 通 量 。 涡 度 相关 及 波 文 比 能 量 平 衡 系统 测 得 的 
通 量 数据 通过 EddyPro ( LI-COR ) 软件 进行 峰值 检 
测 .频率 响应 校正 .声波 温度 转换 为 实际 温度 等 后 台 
处 理 , 并 使 用 平面 拟 合 方法 和 WPL( Webb-Pearman- 
Leuning) 密度 校正 。 使 用 前 均 进 行 坐标 旋转 、 数 据 
插 补 和 求 日 平均 等 处 理 , 以 获得 质量 较 高 的 实测 数 
据 ""”。 利 用 两 类 ET 值 对 SEBS 反 演 的 ET 结果 进 


行 变化 趋势 分 析 及 精度 检验 ,微型 蒸 渗 仪 与 观测 系 
统 如 图 2 所 示 。 


| 
EN 


2 AAR EM 涡 度 相关 (EC) 及 波 文 比 能 量 平衡 ( BREB ) 系统 示意 图 


Fig.2 Schematic diagram of micro-lysimeter ,eddy covariance( EC) and Bowen-ratio-energy-balance( BREB) system 


2 研究 方法 和 模型 


2.1 研究 方法 

2.1.1 SEBS 模型 计算 SEBS 模型 "… ”利用 已 处 
理 完成 的 地 表 参 数 和 气象 数据 ,对 灌区 的 地 表 净 辐 
射 通 量 . 显 热 通 量 及 土壤 热 通 量 进行 估算 。 任 意 时 
刻 的 地 表 能 量 平衡 方程 如 下 : 


R, =G, +hET+H (1) 


AP: R, 为 地 表 获 得 的 净 太 阳 辐 射 通 量 (W， 
m“) ;Cu 为 土壤 热 通 量 (W m°), FERMERE 
的 热 交换 ;NE7 为 潜 热 通 量 (W . m“) ,下 热 面 与 大 
气 之 间 的 水 汽 热 交换 ,其 中 ,入 =2.45 x 10° J kg’, 
为 水 的 汽化 热 ;87 HRE (kg m? s71); HX 
感 热 通 量 , 也 称 为 显 热 通 量 (W . m“) ,下 热 面 与 大 


Tim Tit FE MS FACT 
(1) 净 辐射 (R, ) 公式 为 : 


R,=(1 =a) Roa + Riya — E575 (2) 


式 中 :ea 为 地 表 反 射 率 ;R 为 下 行 太阳 辐射 ;Ri 为 
下 行 长 波 辐射 ;s 为 地 表 比 辆 射 率 ;7 为 传感器 测 得 
的 辐射 地 表 温 度 (C ) ;6 为 斯 蒂 芬 - 波 尔 兹 曼 常 数 
(5.67 x10™), 

(2) 土壤 热 通 量 ( Co ) 公 式 为 : 


Go=R,IT,. +(1-f)(T,-T.)] 


(3) 
式 中 :假定 土壤 热 通 量 和 净 辐 射 通 量 之 间 的 比值 对 
于 完全 植被 覆盖 的 地 区 TT, 为 0.05( 植 被 覆盖 较 好 
KER) ;相对 于 裸 土 设 定 为 工 .为 0.315( 裸 土 区 域 ) ; 
在 这 两 种 极端 条 件 之 间 则 可 以 利用 植被 覆盖 度 ( 矿 ) 
来 进行 插值 计算 ,其 中 大 能够 通过 遥感 数据 来 确 
定 ,通过 下 式 计 算 求 得 。 
_ NDVI - NDVI „i 
e NDVI pa — NDVI „pn 


(3) 显 热 通 量 (五 ) 公 式 为 : 
在 大 气 近 地 层 中 ,根据 大 气 边界 层 相 似 理论 ,有 
以 下 的 关系 : 


(4) 


ot) o 
pC,0,u 

L=- heH 7) 

式 中 :z 为 参考 高 度 (m) su ABE Cm + 87!) su" 为 


摩擦 风速 (ms”) ;do 为 零 平 面 位 移 高 度 (m ) ;zw 
和 zu 分 别 为 动力 学 粗糙 长 度 (m) 和 地 表 热 传输 粗 
BERKE (m); Pa 和 分 别 为 动力 学 和 热力 学 传输 
的 稳定 度 订正 函数 ;9 和 6. 分 别 为 观测 面 和 参考 面 
高 度 的 虚 温 (% ) ;L 为 莫 宁 霍 夫 长 度 (m) ;三 为 感 热 
Wa (W-m*) 让 为 卡尔 曼 常 数 ;p 空气 密度 (g 
m -3) ;C, IZARRI g- C) ;9, 近 地 表 的 位 
HACC) sg 重力 加 速度 (m . s”)。 摩 擦 风速 (u) E 
Pull et (A) ` 莫 宁 霍 夫 长 度 (Z) 可 以 通过 返 代 以 上 3 
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个 求解 方程 得 到 ,其 他 的 变量 通过 气象 观测 信息 结 
合 遥 感 观 测 信 息 求 得 。 

(4) HRE (ET) 计算 公 式 为 : 

假定 在 干燥 地 表 环 境 下 ,由 于 土壤 地 表 水 分 亏 
BE, SPT BEAK ot FR BE ED ,此 时 潜 热 遂 量 
AHE ,这 种 条 件 下 ,根据 地 表 能 量 平 衡 方程 ,此 时 
感 热 通 量 达到 最 大 值 , 即 干燥 条 件 下 的 感 热 通 量 
H: 


Ay, = R, = Co (8 ) 


RP: Ha 为 干燥 地 表 条 件 下 的 感 热 通 量 (W- 
2) ;在 土壤 水 分 充足 时 , 湿 限 条 件 下 的 感 热 通 量 


m`“) 
人 ,是 通过 FAO Penman-Monteith 公式 得 出 : 


et 


Hg = [(,-G) a /[1 +] 


_ 1 fy (2-40 z~d, Zon 
ce A (3) (10) 


3 


pu 


(9) 


Ce ear (Ha 
HP HIRERE A, 定义 为 : 
OAE ©, AEw-AE  H-H, 
a DE ye 7 AE va Ha =H ya tia) 


RRE A FE ASE Ds a BAR H H RE T BY) E 


A, AE, 
= AE 2 we (13) 
R, ~ Co R, 一 Cn 


综合 上 式 ,可 得 到 日 57 计算 公式 为 : 


A 


R, z Go 


ET =8.64 x10’ xA x 
Ap 


(14) 
APET AAA AR BG A 为 蒸发 比 ; p, 为 水 的 密 
度 。 为 了 确定 蒸发 比 的 值 , 需 要 考虑 在 最 干 和 最 湿 
润 两 种 极限 状态 下 的 能 量 平衡 关系 。 

2.1.2 土地 利用 分 类 ”参照 2017 年 哨兵 2 号 全 国 
土地 利用 分 类 图 及 分 类 标准 " ,通过 全 新 的 图 像 分 
类 框架 (ENVI Classification Framework ) 提供 的 土地 
利用 分 类 方法 Softmax ,将 河套 灌区 的 土地 类 型 划分 
为 水 体 耕地、 林地 ,城市 (居民 用 地 、 建 设 用 地 ) 、 章 
HAN Fes 6 类 。Softmax 土地 分 类 方法 使 用 不 同 算 
法 进行 训练 ,充分 利用 对 象 的 光谱 特征 纹理 特征 与 
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空间 特征 等 ,输入 多 波段 组 合 的 栅 格 图 像 和 训练 样 
本 ,在 保证 分 割 对 象 间 的 平均 异 质 性 最 小 和 对 象 内 
像 元 间 同 质 性 最 大 的 基础 上 ,通过 梯度 下 降 法 和 和 迭 
代 法 对 像 元 进行 分 割 与 合并 ,实现 遥感 图 像 不 同 地 
物 类 型 分 类 。 本 研究 利用 损失 曲线 和 混淆 矩阵 对 土 
地 分 类 精度 进行 评定 。 损 失 曲 线 是 该 分 类 方法 对 分 
类 结果 的 初步 判定 。 随 迭代 次 数 增加 ,损失 曲线 逐 
MFA ,变化 趋势 平稳 时 ,可 进行 分 类 评估 和 图 像 
分 类 。 混 涌 和 矩阵 是 一 种 用 来 表示 分 类 精度 的 标准 格 
式 , 是 由 一 组 n xn 的 矩阵 构成 ,表示 分 类 数 ,P; 为 
参考 数据 i 类 被 分 到 分 类 数据 j 类 的 样本 数 ,P; 为 正 
确 分 类 数 ,P;, 为 分 类 得 到 的 第 i 类 总 和 ,P, ,为 分 类 
数据 的 第 7 类 总 和 ,P 为 样本 总 数 。 根 据 混 清 矩阵 
法 对 分 类 结果 精度 评价 ,评价 的 指标 主要 有 生产 者 
精度 (PA) ,用户 精 度 (UA) ,总体 精度 (OA ) K Kap- 
pa RUKIA), 


Pi 
PA =p“ (15) 
UA = Pi (16) 

ao 

EP; 
OAS (17) 

PÈP; -2 (PaPa) 

KIA =—= a (18) 

P 

2.2 模型 验证 方法 


利用 不 同方 式 得 到 的 两 类 ET 值 ,对 SEBS 模型 
反 演 结果 进行 精度 检验 。 一 是 以 临河 试验 地 蒸 渗 仪 
观测 的 E7 数 据 、 分 子 地 涡 度 相关 系统 及 波 文 比 能 
量 平衡 系统 得 到 的 ET 值 作为 实测 检验 数据 ,二 是 
以 PM 综合 方程 计算 结果 作为 理论 检验 数据 ”| 。 
PM 模型 的 计算 公式 如 下 。 


Le a te dh 
NIAAA IV A TFA 


900 
0. 408A (R, -G,) +y T4273% (es =e.) 
A +y(1 +0. 34u, ) 


式 中 :A 为 饱和 水 汽 压 曲 线 斜率 ,kPa + C'S R, 为 
地 表 净 辐射 量 , MJ . m e d ;Gu 为 土壤 热 通 量 ， 
MJ m” -d'y 为 干 湿 表 常数 ,kPa* C7 ;ww 为 2 
m 处 风速 ,ms ies\es 为 饱和 水 汽 压 、 实 际 水 汽 
JE kPa; T 为 平均 温度 ,5 。 


3 结果 与 分 析 


3.1 结果 验证 
3.1.1 SEBS 模型 反 演 结果 验证 根据 实测 数据 对 
ET 反 演 结果 进行 变化 趋势 分 析 及 精度 检验 ( 表 1) 。 
从 表 1 可 知 :研究 区 SEBS 模型 反 演 结果 与 3 种 实 
WM ET 值 变 化 趋势 基本 一 致 ,SEBS 模型 反 演 结果 与 
波 文 比 观测 值 最 接近 ,不 同时 间 的 相对 误差 差异 较 
小 ,相对 误差 平均 值 为 14.93% 。SEBS 模型 的 反 演 
结果 与 蒸 渗 仪 . 涡 度 系统 观测 值 的 相对 误差 平均 值 
分 别 为 16.66% Fil 18. 34% 。 

利用 国际 应 用 较 广 ,公认 的 PM 模型 对 SEBS 
模型 反 演 结果 进行 点 位 精度 验证 ( 表 2)。 由 表 2 可 
以 看 出 :SEBS 模型 与 PM 模型 的 5 个 点 位 ET 值 随 
时 间 变 化 趋势 相 一 致 , 受 不 同 气象 参数 的 影响 ,五 原 
的 ET 估算 结果 与 PM 模型 的 ET 计算 结果 最 接近 ， 
相对 误差 平均 值 为 8. 23% ,9 月 的 相对 误差 最 小 ,为 
1.05% ,5 .6 月 的 相对 误差 接近 且 达 到 最 大 ,相对 误 
差 为 15. 40% 。 临 河 的 ET 估算 结果 与 PM 模型 的 
ET 计算 结果 差异 略 大 ,相对 误差 平均 值 为 20. 94% , 
7 月 两 者 的 ET 结果 较 接 近 , 相 对 误差 为 5. 61% ,5 
月 的 差 值 较 大 ,相对 误差 为 27.61% 。 受 遥感 影像 
质量 与 气象 数据 影响 ,不 同 点 位 的 相对 误差 值 在 时 
间 上 分 布 不 一 致 ,导致 SEBS 模 型 反 演 结果 的 偏离 


ET y = 


(19) 


表 1 AA SEB 模型 计算 值 与 实测 结果 对 比 
Tab.1 Comparison of SEBS modelcalculated values and measured results 
月 份 SEBS 估算 值 波 文 比 观 测 值 相对 误差 涡 度 系统 观测 值 ”相对 误差 ACEO == 相对 误差 
/(mm:d-!) /(mm:d-!) /% /(mm + d7!) 1% /(mm + d7!) 7% 
5 4.52 = = = = 5.58 19. 00 
6 4.72 5.23 9.75 = = 5.64 16.31 
7 5.21 5.98 12.88 6.39 18.47 5.69 8.44 
8 4.27 5.16 17.25 5.36 20.34 5:33 19. 89 
9 3.51 4.38 19.86 3.02 16.23 4.37 19. 68 
相对 误差 平均 值 14.93 18.34 16. 66 
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表 2 不 同 气 象 站 点 SEBS 模型 计算 值 与 PM 计算 值 对 比 
Tab.2 Comparison of SEBS and PM calculated values at different weather stations 


临河 TEH 


杭锦后旗 五 原 


乌拉 特 前 旋 


SEBS PM ”相对 SEBS PM ”相对 SEBS 
月 份 
/(mm+ /(mm* 误差 /(mm: 
d') da) /% dt) da) /% ad) 


/(mm: 误差 /(mm + 


PM ”相对 SEBS PM ”相对 SEBS PM ”相对 
/A(mm*， 误差 /(mm: /(mm: 误差 /(mm* /(m: 误差 
d!) /% d7!) d!) /% d7!) d!) /% 


5 4.022 5.556 27.61 4.788 5.662 15.44 4.641 

6 4.488 5.684 21.04 4.929 6.278 21.49 4.889 

7 5.465 5.790 5.61 5.779 6.978 17.18 6.324 

8 3.819 5.094 25.03 4.643 5.081 8.62 4.431 

9 3.175 4.256 25.40 4.276 4.503 5.04 4.278 
相对 误差 20.94 13:55 


平均 值 


程度 不 同 。PM 模型 计算 的 5 个 气象 站 的 相对 误差 
平均 值 为 14.03% , 均 在 允许 的 误差 范围 之 内 ,SEBS 
模型 估算 整个 河套 灌区 ET 具有 较 高 的 准确 性 和 实 
用 性 。 

3.1.2 土地 利用 分 类 结果 验证 以 哨兵 2 号 土地 
利用 分 类 结果 作为 本 研究 土地 分 类 参照 标准 ,如 图 


5.344 13.15 4.662 5.510 15.39 4.840 5.836 17.07 
5.742 14.86 4.715 5.573 15.40 4.938 6.519 24.25 
5.792 9.19 5.707 5.597 1.97 5.852 6.592 11.23 
5.152 13.99 4.593 4.957 7.34 4.583 5.801 21.00 
4.465 4.19 3.935 3. 894 1.05 4.141 4.512 8.22 

11.08 8.23 16.35 
3a 所 示 。 基 于 Landsat 8 遥感 影像 ,通过 Softmax 分 


类 方法 得 到 灌区 土地 类 型 分 类 结果 (图 3b)。 以 
Softmax 作为 分 类 器 的 土地 利用 分 类 结果 生成 损失 
曲线 与 分 类 混淆 矩阵 ,损失 曲线 如 图 4 所 示 。 通 过 
分 类 混 清 矩阵 进行 精度 评价 ,各 分 类 类 别 的 生产 者 
精度 和 用 户 精 度 如 表 3 所 示 。 


(四 哨兵 ?号 m 耕地 m 草地 mki 加 荒地 


m 林地 m 灌木 林地 m 城市 用 地 


41°N 


01020 40 60 80km 


(b)Landsat 8 mm 水 体 四 城市 用 地 sët 


om 耕地 mm 林地 =e 


41°N 


01020 40 


60 80km 


107°E 108° E 109° E 


107°E 108° E 109° E 


图 3 哨兵 2 号 土地 分 类 与 Landsat8 土地 分 类 


Fig.3 Sentinel No. 2 classification and Landsat8 classification 


20) Saal i eo po ial ee 


West 4) BUA: 0.000 01 
ARR: 100 
最 大 迭代 次 数 : 100 


Loss 


Dl i lin rlia Lt 
0 20 40 60 80 100 
Iterations 


图 4 损失 曲线 
Fig.4 Loss profile 


表 3 土地 利用 分 类 精度 评价 结果 


Tab.3 Land use classification accuracy evaluation results 


城市 Bes 


类 别 水 体 ”耕地 用 地 林地 “荒地 草地 
水 体 2 142 10 516 24 58 198 
耕地 10 3140 762 32 123 254 
城市 用 地 0 0 5 155 64 382 90 
林地 0 0 5 381 0 0 
Fie SH 0 14 1 029 39 7 575 273 
草地 0 3 62 32 54 421 
分 类 总 数 2152 3167 7529 572 8192 1236 
生产 者 精度 
A 99.54 98.83 78.47 75.09 92.47 76.50 
用 户 精度 
(UA/% ) 67.19 70.80 90.42 66.87 84.54 73.61 


总 体 精度 (OA) =80.63% Kappa 系数 (KIA) =0.75 
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通过 损失 曲线 调整 分 类 学 习 参 数 ,得 到 最 优 分 
类 结果 。 损 失 曲 线 逐 渐 趋 于 零 , 量 变化 趋势 平稳 时 ， 
可 用 于 本 次 研究 分 类 评估 和 图 像 分 类 。 基 于 Land- 
sat8 影像 土地 利用 分 类 结果 与 哨兵 2 号 分 类 结果 的 
空间 分 布 基本 一 致 ,Softmax 作为 分 类 器 的 总 体 分 类 
精度 为 80. 63% ,Kappa 系数 为 0. 75 ,分 类 可 信和 度 较 
高 。 治 区 内 林地 只 分 布 在 道路 两 侧 ,导致 林地 的 类 
别 总 和 数 较 小 且 分 类 精度 较 低 。 城 市 用 地 中 夹杂 小 
面积 其 他 土地 类 型 ,导致 分 类 精度 降低 。 
3.2 SEBS 模型 估算 ET 结果 分 析 
3.2.1 A ETH MRE 灌区 生长 季 ETR 
结果 如 图 5 所 示 ,E7T 值 从 5 月 开始 逐渐 上 升 ,7 月 
达到 最 大 值 ,7 一 9 月 逐渐 下 降 。 总 体 呈 单 峰 形 分 
布 ,峰值 出 现在 7 月 。 整 个 河套 灌区 7 HH ETH 
5.21 mm-d',6 月 与 5 月 的 日 87 分 别 为 4.72 
mm"d 4.52 mm .d ,8 月 与 9 月 的 日 £7 分 别 
为 4.27 mm + d7' 3.51 mm: d7', 
3.2.2 日 ET 空间 差异 特征 结合 河套 灌区 ET 时 
空 分 布 图 (图 5) 与 土地 利用 分 类 图 (图 3) ,不 同 土 
地 类 型 ET 随时 间 也 会 发 生 相 应 变化 。 灌 区 生长 季 
内 ,水 体 湖区 水 量 充沛 ,因此 5 一 9 月 的 水 体 湖区 ET 
值 较 大 ,主要 集中 在 乌 梁 素 海 .金马 湖 等 。 其 次 耕地 
ET 主要 影响 因素 是 人 为 干预 ,灌区 定期 灌溉 导致 种 
植 作 物 的 耕地 区 域 有 充沛 的 水 分 供 植物 蒸腾 与 土壤 
蒸发 。 草 地 略 大 于 荒地 的 E7 ,与 荒地 相 比 ,草地 有 
来 自 周 边 湖区 或 耕地 区 灌溉 后 少量 用 水 、 季 节 性 降 


水 所 提供 的 水 源 补给 。 荒 地 含 盐 量 较 大 且 储 水 能 
低 ,不 适宜 作物 生长 ,导致 长 期 荒 置 ,ET 变化 较 小 ， 
处 于 极 低 值 ,灌区 的 西南 方向 邻近 三口 县 的 区 域 万 
为 明显 。 城 市 区 域内 一 般 夹 杂 小 面积 林地 ,水域 . 草 
地 等 用 地 类 型 ,不同 月份 的 城市 用 地 ET 值 均 处 于 
中 等 水 平 ,随时 间 变 化 不 明显 。ET 高 低 区 域 间 具 有 
明显 的 分 界 ,与 不 同 土地 类 型 的 边界 较 一 致 ,表明 不 
同 土地 类 型 ET 差异 明显 ,所 需 水 资源 .土地 利用 方 
式 也 不 相同 ,可 以 用 来 指导 灌区 内 水 资源 管理 分 配 
与 土地 利用 规划 。 
3.3 不同 土地 类 型 ET 差异 分 析 

河套 灌区 受 人 类 活动 强烈 影响 ,ET 也 受 很 多 
素 影响 ,不 同 的 土地 类 型 ET 也 会 存在 明显 差异 。 
为 进一步 分 析 不 同 地 物 覆 盖 类 型 下 的 ET 变化 规 
IE ,将 灌区 划分 为 水 体 、 耕 地 、 林 地 城市 (居民 用 
地 建设 用 地 ) 草地 和 荒地 6 大 类 ,不 同 土 地 类 型 
的 面积 比例 如 表 4 所 示 , 可 以 看 出 灌区 主要 由 耕地 、 
草地 .荒地 组 成 ,占据 灌区 总 面积 的 93% ,水 体 和 林 
地 的 面积 最 小 , 仅 占 3% ,城市 用 地 占 4% 。 

从 图 6 可 以 看 出 ,ET7 值 从 5 月 作物 生长 开始 到 
7 月 逐渐 升 高 ,7 一 9 月 逐渐 降低 ,灌区 6 类 土地 类 型 
的 ET 变化 趋势 基本 一 致 。 水 体 与 其 他 土地 类 型 的 
下 垫 面 不 同 , 有 充足 的 水 分 蒸发 ,因此 以 乌 梁 素 海 为 
典型 代表 的 水 体 区 域 ET 最 大 。 林 地 ET 仅 次 于 水 
体 ,土壤 水 分 长 期 较 大 ,树木 不 断 进行 蒸腾 蒸发 。 耕 
地 ET 的 变化 最 为 明显 ,6 一 7 月 耕地 的 87 急 剧 增 


(a) 2017-05 (b) 2017-06 
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Fig.5 Daily evapotranspiration results for the study area on a monthly basis 
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表 4 不 同 土地 利用 类 型 面积 比例 


Tab.4 Percentage areas of different land use types 


土地 利用 类 型 ”水体 ”耕地 城市 用 地 林地 荒地 ”草地 
面积 比例 /% 1.98 74.78 3.94 0.73 12.67 5.90 
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图 6 不 同 土地 利用 类 型 蒸 散 对 比 


Fig.6 Evapotranspiration comparisons for different 


land use types 


大 ,是 人 为 因素 对 其 进行 了 大 面积 .水量 充足 的 灌溉 
所 致 。 草 地 ,城市 用 地 sie a = AY) ET 变化 趋势 基 
本 一 致 ,草地 与 荒地 的 ET 值 基本 相同 ,草地 ET 仅 
略 高 于 荒地 0. 1 mm“'d ,两 者 无 外 界 直接 干扰 ， 
ET 主要 取决 于 降雨 .土壤 水 分 及 地 下 水 等 。 城 市 用 
地 夹杂 不 同 的 土地 类 型 ,57 值 仅 受 降水 影响 ,与 草 
地 ,基地 的 £7 差异 较 小 。 


4 结论 


本 研究 基于 地 表 能 量 平衡 系统 建立 的 SEBS 模 
型 ,以 气象 数据 及 Landsat 8 遥感 影像 为 源 数 据 ,得 
到 内 蒙古 河套 灌区 2017 年 5 一 9 月 的 日 ET 估算 结 
果 。 将 藻 渗 仪 观测 数据 、 涡 度 相关 及 波 文 比 实测 数 
据 作为 实测 值 ,PM 模型 计算 结果 作为 理论 值 ,对 ET 
反 演 结果 进行 变化 趋势 分 析 及 精度 检验 。 土 地 利用 
分 类 方法 Softmax 将 研究 区 分 为 6 类 ,结合 模型 估算 
结果 ,对 不 同 土地 类 型 ET 值 进 行 差异 分 析 。 主 要 
结论 如 下 : 

(1) 实测 数据 和 PM 计算 值 对 SEBS 估算 结 
精度 验证 及 变化 进行 趋势 分 析 , 估 算 结 果 的 变化 趋 
势 与 4 种 验证 数据 变化 趋势 相 一 致 ,相对 误差 在 允 
许 范 围 内 ,结果 对 于 研究 河套 灌区 的 不 同 土地 类 型 
及 区 域 57 具有 较 高 应 用 价值 。 

(2) KEKEN H ET AIREZ, ETW 
5 月 开始 逐渐 增 大 ,到 7 月 达到 最 大 值 ,日 ET 值 达 
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5.21 mm . d, 7—9 月 逐渐 减 小 ,9 AA ET 值 最 
小 ,为 3.51 mm + do! ET 值 大 小 随时 间 变 化 呈正 

(3) 结合 灌区 土地 类 型 分 布 与 57 值 空间 分 布 
得 知 ,不同 土地 利用 类 型 的 日 ET 大 小 依次 为 :水 
体 > 林地 > 耕地 > 草地 > 城市 > 荒地 ,ET 值 依次 
为 :5.007 mm + d 4.904 mm - d7' 4.742 mm - 
d 4.694 mm + d™' 4.670 mm + d7' 4.643 mm - 
d ,ET 差异 值 分 界 与 灌区 土地 类 型 的 分 布 高 度 一 
致 。 

(4) 生长 季 的 不 同时 间 , 耕 地 的 57 变化 趋势 
与 作物 生长 情况 相 一 致 ,水 体 、 林 地 的 ET 长 期 处 于 
较 高 值 , 草 地 ,城市 用 地 与 荒地 的 ET 较 小 , 且 三 者 
AY ET 值 相 接近 。 
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Estimation and variability of evapotranspiration for different land types 


during the growing season in the Hetao Irrigation District 


WANG Yan-xin', LI Rui-ping', LI Xia-zi 
(1. Water Conservancy and Civil Engineeling College ,Inner Mongolia Agricultural University , Hohhot 010018 , Inner Mongolia , China ; 
2. Insititue of Meteorological Research of Inner Mongolia , Hohhot 010051 ,Jnner Mongolia , China) 


Abstract: To reveal the differences in evapotranspiration (ET) across different land types in the irrigation dis- 
tricts ,the single-day ET during the 2017 growing season in the Hetao Irrigation District was estimated based on the 
surface energy balance system (SEBS) model, meteorological data ,remote sensing images, and land use classifica- 
tion results for the irrigation districts. The SEBS model accuracy was verified using the measured ET data and the 
calculated values obtained using the United Nations Food and Agriculture Organization Penman-Monteith ( PM) 
model. The objective of this study was to analyze the temporal and spatial changes in ET in the irrigation area and 
the differences in ET across different land types. The results showed that (1) the average relative errors obtained 
based on the PM model, eddy covariance , Bowen ratio-energy balance system, and micro-lysimeter observation data 
were 14.03% ,18. 34% ,14.93% ,and 16.66% , respectively. The SEBS model exhibited a high application poten- 
tial for studying the ET across the Hetao Irrigation District. (2) During the growing season, the daily ET values in 
the irrigation area exhibited the following relation; July > June > May > August > September. The daily ET in 
the irrigation area had a unimodal distribution and reached the maximum value of 5.21 mm .+ d~' in July. The daily 
average ET was the lowest in September at 3.51 mm + d~'. (3) The daily ET values for different land types exhibi- 
ted the following relation: water > forest > cropland > grassland > urban > bare land. The ET variability in the 
irrigation area was consistent with the land type. The overall classification accuracy and kappa coefficient for land 
use classification were 80.63% and 0.75 ,respectively. (4) The ET value for water was the higher during different 
periods of the growing season; the daily ET reached the maximum value of 6.034 mm + d~' in July, whereas the 
daily average ET was the lowest in September at 4.177 mm - d~'. The ET value for forest within the irrigation area 
was relatively large, with a daily ET of 5.977 mm + d7' in July and 4. 147 mm - d™' in September. The ET value 
for cropland changed with the growth of the crops,reaching 5.851 mm .+ d~' during the vigorous growing period and 
3.952 mm - d during the harvesting period. 

Key words: growing season; evapotranspiration; SEBS model; Softmax classification method; land types; Hetao 


Irrigation District 


